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INTRODUCTION 


The  formation  of  the  nocturnal  inversion  was  probably  first  examined 
formally  by  Taylor  [1].  His  investigation  was  based  upon  the  eddy 
transfer  of  heat  with  the  eddy  conductivity  independent  of  height, 
which  can  be  stated  as 


3T  _ ^ 
3t  ” ' 


32T 
H 32 


(1) 


where  T is  temperature,  t time,  and  z height  above  the  surface.  For  the 
case  of  a nocturnal  inversion,  where  the  rate  of  temperature  change  at 
the  surface  can  be  considered  to  be  fairly  uniform,  Taylor  suggested 
that  the  temperature  diminishes  n-degrees  per  unit  time  so  that  at  time 
t the  temperature,  T,  is  T^  - nT  where  T^  is  the  temperature  at  t^. 

At  time  t,  the  temperature  at  z will  be 


where  6 is  a constant  lapse  rate,  and  p the  kinematic  viscosity.  Taylor 
found  that  the  term  multiplying  nt  was  unity  at  the  surface  and  0.1  at 
z(4K^t)"**  = 0.8  and  suggested  that  no  effect  existed  beyond  z(4Kj^t)"^  = 1, 
so  that  surface  changes  of  temperature  over  time  t extend  to  a height 
defined  by 


zf  = 4K„t 


(3) 


where  z^  is  the  top  of  the  inversion. 

Sutton  [2],  using  a multilayer  approach,  but  similar  reasoning,  verifies 
the  approximation  of  Eq.  (3)  for  the  growth  of  nocturnal  inversions  with 
time  where  Ku  is  a maximum  at  some  height  h above  the  surface. 
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Estimates  of  the  depth  of  nocturnal  inversions  are  of  some  importance  to 
transport  and  dispersion  of  smoke/obscuration  materials  released  into  the 
atmosphere  during  stable  conditions.  Crosswind  integrated  concentrations, 
source  strengths,  and  downwind  diffusion  of  materials  released  into  the 
atmosphere  are  highly  dependent  upon  stability  and  inversion  depths.  In 
turn,  these  parameters  determine  munition  expenditures  and  obscuration. 

BASIC  CONSIDERATIONS  FOR  THE  SURFACE  BOUNDARY  LAYER 
The  lower  portion  of  the  atmosphere  may  be  considered  to  consist  of  two 
layers,  the  surface  boundary  layer  < z < h and  the  planetary  layer 

proper  extending  to  the  gradient  wind  level  h < z < Z,  where  z^  is  the 
roughness  length,  h the  top  of  the  surface  layer,  and  z the  gradient 
wind  level.  It  may  also  be  assumed  that  the  exchange  coefficients  for 
heat  and  momentum  reach  a maximum  at  h,  decreasing  again  to  a residual 
value  at  z. 

In  a thermally  stratified  stable  regime,  vertical  fluxes  of  heat  and 
momentum  can  be  considered  to  be  functions  of  mechanical  turbulence  alone, 
which  suggests  that  where  is  the  eddy  viscosity,  with  the 

dimensionless  parameters  of  the  Obukhov  [3]  similarity  theory  given  by 


and  Rf  = 


where  L is  the  Obukhov  scaling  length,  Ri  and  R^  the  gradient  and  flux 
Richardson  numbers,  respectively,  with  ^ a dimensionless  wind  shear 
and  a dimensionless  lapse  rate.  Hansen  [4]  has  demonstrated  that  for 
stable  flow  the  wind  profile  may  be  written  as  i 
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where  V Is  the  mean  horizontal  windspeed,  u*  a friction  velocity,  and  k 
Karman's  constant.  Furthermore,  Hansen  [4]  also  found  that 


^ = R1  + 15  Ri2  ( 

^ = 1+  B^=1+15r1  ( 

where  s Is  the  average  over  a height  Interval  of  a variable  given  by 

® ■ "Ri  7/1'  ^ 


The  mean  reduces  to 


^ “ rr-H-RT 


If  the  wind  profile  In  finite  difference  form  Is  written  for  stable  flow 


AV 

imrnr  '■ 


where  z ■ h « L R1  0jy|,  then 

If  aV  Aln'^z,  the  critical  gradient  at  the  geometric  mean  height  h,  is 
assumed  to  occur  over  the  layer  Ain  z ® In  e » 1,  then  RIaV  = , and 

h » mL  = ml  (11 

where  m Is  the  profile  slope  k/u*.  If  h/L  = R1  then  from  Eq.  (11) 

R1(h)  • * (1! 


and  from  Fq.  (5) 


h = L (R1(h)  + 15  Ri(h)2). 


(13) 


The  height  h is  taken  to  be  the  depth  of  the  surface  boundary  layer  in  a 
stable  flow  regime  and  is  the  height  where  or  maximum  is  evaluated 
for  inversion  depth  estimates  using  Eq.  (3).  The  eddy  viscosity  at  h is 
given  by 


ku*z 


(14) 


Lumley  and  Panofsky  [5]  suggest  that  the  depth  of  the  surface  layer  can 
be  estimated  from 


h = 20  Tq  (15) 

where  is  the  surface  stress  and  assumed  to  be  the  equivalent  of  the 
Reynolds  shearing  stress,  t = puj,  where  p is  density.  Calculations 
based  upon  stable  regime  data  extracted  from  Lettau  and  Davidson  [6]  and 
Barad  [7],  as  sutimarized  in  Figure  1,  indicate  that  Eq.  (15)  provides 
estimates  of  h comparable  with  Eq.  (13). 

ESTIMATES  OF  INVERSION  HEIGHTS 

The  scheme  for  estimating  inversion  heights  was  evaluated  by  using  experi- 
mental data  extracted  from  Lettau  and  Davidson  [6];  Barad  [7];  Bowne, 
Entrekin,  and  Smith  [8];  and  Stenmark  and  Drury  [9].  Micrometeorological 
profile  data  from  Lettau  and  Davidson,  Stenmark  and  Drury,  and  Barad  were 
utilized  to  calculate  h,  and  Zj  for  the  thermally  stratified  stable 
regime.  Rawinsonde,  aircraft  soundings  and  wiresonde  data  summaries  from 
all  four  data  samples  were  used  to  obtain  indicated  inversion  heights. 
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According  to  Milly  [10],  inversion  conditions  can  be  assumed  to  exist 
from  approximately  60  minutes  before  sundown  to  about  60  minutes  after 
sunrise.  For  those  data  samples  where  sunrise  and  sunset  were  not  listed, 
these  times  were  estimated  from  information  in  The  American  Ephemeris 
and  Nautical  Almanac,  issued  by  the  Nautical  Almanac  Office,  US  Navy 
Observatory,  Washington,  DC. 

The  calculated  and  observed  values  for  Zj  from  the  217  micrometeorological 
profiles  and  the  305  soundings  used  in  this  study  were  averaged  by  using 
two  methods:  (1)  hourly,  and  (2)  by  means  of  a geometric  progression 
where  the  overlapping  time  intervals  serve  to  filter  and  smooth  the  data. 
The  hourly  averages  of  Zj  (observed)  and  Zj  (calculated)  are  shown  in 
Figure  2 and  the  overlapping  time  averaged  analysis  in  Figure  3.  Note 
that  the  calculated  values  of  Zj  overpredict  inversion  heights  soon  after 
sunset  and  underpredict  in  the  early  morning  hours.  No  significant  dif- 
ference is  apparent  between  the  two  averaging  methods. 

To  simplify  the  prediction  of  inversion  growth  and  height  with  time,  a 
semiempirical  formula  based  upon  Eq.  (3)  where  the  value  of  Kjj|(h) 
estimated  from  0^  and  V was  developed  and  takes  the  form 

Zj  = 7.750-^  (V  z t)^  (16) 

where  z is  2 m,  7 is  windspeed  in  m sec"^  and  t = minutes.  Figure  4 
shows  values  of  z^  calculated  by  using  Eq.  (16)  with  respect  to  the 
observed  average  inversion  heights  with  respect  to  time. 
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CONCLUSIONS 

The  results  shown  in  Figures  2 through  4 are  in  good  agreement  with  a 
study  reported  on  by  the  US  Weather  Bureau  [11]  as  summarized  by  Wanta  [12]. 
It  should  be  pointed  out  that  the  conditions  indicated  in  Figures  2 
through  4 are  averages  based  upon  four  distinctly  separate  data  samples 
observed  in  widely  separated  geographical  locations  over  dissimilar  ter- 
rain. The  variations  in  observed  or  calculated  inversion  heights,  as  indi- 
cated by  the  one  standard  deviation  bars  about  the  meaned  values  of  the 
figures,  represent  the  variations  that  could  occur  in  inversion  heights, 
even  in  a single  nocturnal  period.  This  can  be  verified  by  considering 
the  remotely  sensed  data  of  Hall  [13]  and  the  summary  of  the  status  of 
remote  sensing  prepared  by  Little  [5].  The  fluctuations  in  the  inversion 
height  with  time  can  be  easily  recognized  in  Little's  Figures  10  and  11. 

One  difficulty  encountered  while  operating  upon  the  profile  data  using 
Eqs.  (3)  and  (16)  was  the  low  windspeed  cases.  Speeds  less  than  2 to  3 
m sec”^  at  2 meters  height  tended  to  underestimate  inversion  heights 
badly,  especially  6 to  8 hours  into  the  nocturnal  regime.  Thus,  the  use 
of  the  prediction  estimators  is  not  recomnended  during  periods  of  low 
windspeeds.  Generally  though,  Eqs.  (3)  and  (16)  do  a reasonably  accurate 
job  of  estimating  inversion  heights  as  long  as  the  atmosphere  is  stable 
and  the  2-meter  windspeed  is  greater  than  3 m sec"^. 

The  results  presented  are  applicable  to  US  Army  problems  associated  with 
the  dispersion  and  transport  of  smoke/aerosols  utilized  for  screening  and 
obscuration. 
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